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Abstract. We revisit the problem of finding the probability distribution of a 
fermionic number of one-dimensional spinless free fermions on a segment of a given 
length. The generating function for this probability distribution can be expressed as a 
determinant of a Toeplitz matrix. We use the recently proven generalized Fisher- 
Hartwig conjecture on the asymptotic behavior of such determinants to find the 
generating function for the full counting statistics of fermions on a line segment. Unlike 
the method of bosonization, the Fisher-Hartwig formula correctly takes into account 
the discreteness of charge. Furthermore, we check numerically the precision of the 
generalized Fisher-Hartwig formula, find that it has a higher precision than rigorously 
proven so far, and conjecture the form of the next-order correction to the existing 
formula. 
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1. Introduction 

The problem of finding the fermion-number fluctuation on a segment for free one- 
dimensional spinless fermions is a classical textbook problem. It is related to a variety of 
topics in the theory of one-dimensional quantum systems and in mathematical physics, 
such as quantum spin chains [HE], bosonization (3], full counting statistics (FCS) [HE], 
random-matrix theory [BJ, theory of Toeplitz determinants [7J [S], etc. Remarkably, 
some helpful results related to this simple problem have been obtained only recently: 
a generalized version of the Fisher-Hartwig conjecture on Toeplitz determinants has 
been proven p]. These mathematical results were further used in a non-equilibrium 
bosonization approach in Refs. [Till ITT]. 

In the present paper, we use the generalized Fisher-Hartwig formula to calculate the 
generating function for FCS of free fermions on a one-dimensional line segment. Unlike 
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the bosonization approach, the generalized Fisher-Hartwig formula treats correctly 
the discreteness of particles. As a result, we obtain a quantitative description of the 
development of the singularity in the FCS generating function in the limit of a large 
segment length (a FCS phase transition, in the definition of Ref. |T2]). 

Furthermore, we investigate numerically the precision of the existing generalized 
Fisher-Hartwig formula and come to the conclusion that it is more precise than formally 
proven. In addition, our numerical results allow us to extract the next-order correction 
to the existing formula and to conjecture its analytic form. 

The paper is organized as follows. In Section [2] we define the problem of counting 
free fermions on a one-dimensional line segment. In Section [3j we briefly review the 
existing analytical approaches to the problem and apply the generalized Fisher-Hartwig 
formula to obtain a good approximation for the FCS generating function. In Section |4| 
we analyze the precision of the obtained expression numerically and conjecture its 
improvement by including a higher-order correction. Finally, in Section [5] we summarize 
our findings. Some technical details are presented in the Appendices. 



2. Formulation of the problem: counting fermions on ID line segment 

Consider free spinless fermions in one dimension (either on a line or on a one-dimensional 
lattice) at zero temperature. We study the fluctuations of the fermion number on 
a segment of a large length L. One can easily find (using either Wick theorem or 
bosonization) that 

((Q 2 )) = {{Q-(Q)) 2 )~-ML/lo) (1) 



7T 



in the limit of a large segment length L ^> Iq. Here Q = dxc^c is the operator of 
the total number of fermions (total charge) on the segment (c^ and c are the fermionic 
creation and annihilation operators). The averaging is performed over the ground state 
of free fermions characterized by the Fermi wave vector kp- The ultraviolet cutoff Iq is 
given by l ~ k F x on the line and by / ~ (sin/cp)" 1 on the lattice. The average charge 
on the segment is given exactly by 

(Q) = Lkp/ 7r , (2) 

(so that the average density of fermions is kpjix). 

A more complicated problem is to find all moments of charge (Q k ) or, equivalently, 
the full distribution of probabilities of having a given charge q on a segment. Instead 
of calculating all moments separately, it is convenient to introduce the characteristic 
function 

oo 

x(A) = (e UQ ) = ^P g e^ (3) 

9=0 

also known as the full-counting-statistics (FCS) generating function. The charge 
cumulants may be expressed as its logarithmic derivatives: 

((Q k )} = (-id x ) k logx(\) ■ (4) 

A — 
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It is obvious from the definition of ^ that (i) x(A) is normalized so that x(0) = 1, 
(ii) x(— A) = x*(A), (hi) it is periodic x(A + 27r) = x(A). The latter property follows from 
the fact that the charge operator is integer- valued (all charges are integer numbers). Our 
goal is to find a good approximation for x(A) in the limit of large L. 



3. From bosonization to generalized Fisher— Hartwig conjecture: early 
approaches and new results 

The simplest approach to calculating x(A) is bosonization, which is equivalent to 
assuming that the density fluctuations are Gaussian. This assumption leads immediately 
to the following result 

\nx(X)^^X(Q}-^((Q 2 }} (5) 



with the average charge and the variance given by ^ and ([I]), respectively (see Appendix 
[A] for a brief review of the bosonization calculation). 

One can notice two obvious drawbacks of this approximation. First, the 
bosonization method does not allow to calculate the numerical coefficient for the 
ultraviolet cut-off Iq. This problem can be partly solved by a direct calculation of 
the second moment in the original fermionic representation using the Wick theorem. 
Such a calculation (see Appendix D) reproduces the result ([I]) with 



/- 1 = 2e 7B+1 sinfc F (6) 

in the lattice problem and Iq 1 = 2e lE+l kp in the continuous case (note that the formulas 
for the continuous case may be obtained from the results on the lattice by taking the 
limit kp — > while keeping the product kpL fixed). Here je = 0.57721... is the 
Euler-Mascheroni constant. 

Nevertheless, even after fixing the numerical coefficient in Iq, the bosonization result 
for x(A) is only precise up to a numerical A-dependent coefficient of order one. These 
A-dependent corrections may be interpreted as cumulants of order higher than two, 
which are not captured by bosonization. If we include those corrections, we arrive at 
the following expansion: 

In x(X) = iX^L - In y + F (X) + o(l) (7) 

as L — > oo. The function F (X) may, in principle, be reconstructed from the cumulants 
((Q k )) calculated with the help of the Wick theorem. Such calculations appear to 
be very tedious. Fortunately, the exact expression for F (X) can be obtained using 
the determinant representation of (|3j) and the Fisher-Hartwig formula for Toeplitz 



determinants (see Appendix B for details and references) 

G ' 1 + s) G ( 1 -s 



F {X) = 2 In 



+ ^(7^ + 1), (8) 
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where G(z) is the Barnes G-function [13] . Remarkably, F (X) does not depend on kjr, 
even in the lattice version. [|] 

The second deficiency of the bosonization approximation is that x(A) does 
not obey the periodicity property. This is due to the nature of the bosonization 
approximation which treats the fermionic density as a continuous field and thus ignores 
the discreteness of the fermionic charge. This problem was addressed in Ref. [7J, 
where a phenomenological formula was proposed to restore the periodicity (taking into 
account umklapp processes, in the bosonization language). Aristov's ansatz in Ref. [7] 
corresponds to neglecting the correction F Q (X) with its nontrivial A dependence. 

In the present work, we show that a rigorous formula restoring the periodicity of 
x(A) can be obtained by an application of the recently proven generalized Fisher-Hartwig 
conjecture [5]. It turns out that the correct recipe for the periodic extension of Eq. ^ 
is simply to add two such expressions with shifted values of A in the vicinity of the 
"switching points" A = (2k + 1)tt: 

X(A) ~ Xo(A - 2ku) + Xo (A - 2[k + 1]tt) , (9) 

where lnxo(A) is given by the right-hand side of Eq. ([7]). The two terms are of the same 
order of magnitude at the "switching point" A = (2k + l)n, but one of them becomes 
sub leading (in L) away from this point. For the same reason, only two shifted values of 
A are of relevance at each "switching point": shifts by higher multiples of 2tt produce 
terms decaying as higher powers of L and therefore may be neglected. The details of the 



application of the Fisher-Hartwig formula to our problem are presented in Appendix 

El 



4. Improving the generalized Fisher-Hartwig formula: numerical analysis 

In the proof of the generalized Fisher-Hartwig conjecture in Ref. [9], its precision is 
only estimated as a relative o(l) as L — > oo. This implies that, in the approximation 
i^m, both terms are within the proven precision only exactly at the switching point (and 
therefore, at the switching point, the expansion (J7J) fails). Away from the switching 
point, the subleading term is already beyond the rigorously proven precision (and thus 
the expansion ^ is the best proven estimate). 

In order to rectify this situation, we perform a numerical analysis of the generalized 
Fischer-Hartwig formula ^ based on exact evaluation of Toeplitz determinants of sizes 
up to 5000. We claim that the formula ^ has a higher precision than that rigorously 
proven: the subleading of the two terms in Eq. (|9j) provides the main correction to the 
asymptotic behavior ([7]). Moreover, the next-order correction may be captured with 
the use of the following conjectured formulas (to simplify notation, we specify to the 
interval A G [0,27r]): 

x(A) = Xi(A) + Xi(A-2tt)+ £ , (10) 

| In the continuum limit this expression is identical to Eqs. (4.6) and (4.19) of Ref. [2] obtained from 
the asymptotics of a Fredholm determinant. 
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where 



Xi(A) = exp 



i\—L - — ln- 
ii 2ir z I 



o 



+ F (X) + F 1 (k F ,X)L 



-i 



and the higher-order correction e is of the relative order L 2 . More precisely, it can be 
estimated as 



IXi(A)| + |xi(A-2tt)|J-0(L- 2 ). (12) 
We confirm this conjecture numerically by extracting the coefficient Fi(k F , A) and 



verifying the estimate (12) for the deviation from the fit. As a result of the fit, we find 



that the coefficient Fi(k F , A) is purely imaginary and an odd function of A. Moreover, we 
observe that, to a very high precision, Fi(kF, A) can be described by a simple formula, 

Fi(k F , A) = - - cotk F (13) 



7T 



sec 



Appendix C for details of the numerical calculations). We therefore conjecture that 



this formula is analytically exact for our expansion ( 10 )— ( 12 ) 



To support this conjecture, we have calculated the second and the third cumulants 



to the order L 1 (see Appendix D). While ((Q 2 )) does not produce any contribution at 
the order L -1 , the third cumulant, to the leading order, is given by 

«Q 3 »~^c U F , (14) 



2tt 3 L 



which is consistent with Eq. (13). Furthermore, if there are corrections to Eq. (13), then 
they are of the order higher than A 3 . Since the relative deviation of the numerically 



extracted Fi(k F , A) from the formula (13) does not exceed 10~ 4 to 1CT 5 in the range of 
A = 0.3 ... 1.3 (error bars are higher for A close to zero and to tt), we conclude that the 
numerical coefficients at such corrections, if nonzero, would be below 10 -4 (at least, for 
the values of k F tested). It seems therefore likely that such corrections vanish identically 



and our suggested expression (13) is exact. It would be interesting to verify or disprove 



this conjecture by the methods of Ref. [9], and we leave this for future studies. 



5. Summary and discussion 

The results of this work are twofold. First, we apply the recently proven generalized 
Fisher-Hartwig conjecture to describe the asymptotic behavior of the FCS problem of 
counting free fermions on a line segment. Second, we analyze numerically the precision 
of the generalized Fisher-Hartwig formula and find that, at least in our case, it is more 
accurate than rigorously proven. Furthermore, we numerically extract the next-order 
correction to the generalized Fisher-Hartwig formula and conjecture a simple analytical 
expression for it. 

We hope that our results will be useful in several respects. On the mathematical 
side, we hope that they will stimulate further studies of the asymptotic behavior of 



Toeplitz determinants. In particular, it seems plausible that our formulas (10) and (11) 
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are the beginning of a more general exact expansion in powers of see Eqs. (E.l) 
and (E.2) of Appendix E Another interesting question would be extending the Fisher- 
Hartwig prescription ^ to block- Toeplitz determinants [151 EE E] • Such an extension 
would be relevant for various FCS problems with noninteracting fermions, where the 
generating function is given by the Levitov-Lesovik determinant formula (H [5] related 
to block- Toeplitz matrices. 

On the physical side, the example we consider in this work is probably the simplest 
quantum example of the "nonanalytic" FCS phase, as defined in Ref. |12j . where 
the thermodynamic limit (L — > oo) can be studied in detail. Remarkably, the "sum 
prescription" (|9]), which describes the development of the singularity in the generating 
function, is the same as in FCS phase transitions in classical Markov processes (see, 
e.g., the discussion of the "weather model" in Ref. [12])- It would thus be interesting 
to explore the extent of universality of the formula ([9]), in particular in the case of 
interacting quantum models. 
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Appendix A. Bosonization calculation of FCS for free fermions 

Here we briefly derive ^ (together with (|2| and ([I])) using the bosonization technique 
(see, e.g., Ref. [3]). In the simplest version of the bosonization approach, the density of 
one-dimensional fermions is given by 

p{x) « p - -d x <p{x), (A.l) 

7T 

where <fi{x) is a free bosonic field with the correlation function 

([0(x)-0(O)] 2 )^lnM (A.2) 



with l G being an ultraviolet cutoff (of the order of lattice constant). The expression (A.2) 
is good at \x\ ^> Iq. The number of fermions on a segment of a length L is obtained 
from (A.l ) as 

Q = J dxp{x) =p o £-^(0(£)-0(O)). (A.3) 
Upon substituting this bosonized form in ([3]), we find 

= e iXpoL (L/l y X2/{2n2) , (A.4) 
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which gives us 

This calculation misses the most important property of its periodicity in 



A. This is due to the approximation (A.l), which neglects terms oscillating with 



the wavevector 2kp and its multiples (and thus carries no information about the 
particle discreteness). The effect of such terms was taken into account in Ref. [7J at 
a phenomenological level, which resulted in a qualitatively (but not quantitatively) 
correct approximation. The ansatz of Ref. [7J can be also reproduced in a different 
way: suppose we calculate the probabilities of different particle numbers as the Fourier 



transform of the bosonization result (A. 4) and then restrict the particle number to be 



integer. This prescription produces the particle-number probabilities 

P q cxe 2 <« 2 » , (A.5) 



which result in the generating function 

+oo +oo 2 

X (A) = J2 P ^ Xq ^ E e^- 2 ^-^^^ , (A.6) 

q=— oo j = — oo 

where we used Poisson's summation formula (the overall coefficient being determined 
from the normalization condition x(A = 0) = 1). Remarkably, this ansatz coincides with 
the Aristov's conjecture in Ref. [TJ. It is only qualitatively correct (the errors being of 
relative order 0(L°)), as can be seen from the exact results based on the Fisher-Hartwig 
formula. 



Appendix B. Fisher-Hartwig formula for FCS of free fermions 

In this Appendix, we review the Toeplitz-determinant approach to the problem and 
apply the generalized Fisher-Hartwig formula. 



Appendix B.l. Toeplitz-determinant representation for x(A) 

We start by reproducing the Toeplitz-determinant representation of x(A) following pp. 
We consider lattice fermions hopping on an infinite one-dimensional lattice with sites 
labeled by integer numbers. The single-particle states are parameterized by the wave 
vector k from the Brillouin zone, k G [— 7T, n]. We assume that the ground state consists 
of the filled states with wavevectors — kp < k < kp, where kp is the Fermi wavevector. 
We do not specify Hamiltonian here, as we are interested in the static problem, where 
the correlator (|3]) is determined only by the ground state of the system. The one-particle 
Green's function in the coordinate space is given by 

-g(k)e-^\ (B.l) 

where the Green's function in momentum space is the step function equal to 1 for 
| A; | < k F and otherwise: 

g(k) = Q(k F - \k\) (B.2) 
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Explicitly combining (B.l) and (B.2), we obtain 



9 



sm(k F (i-j)) 



1-3 



(B-3) 



with g = kp/ir. The operator of the number of fermions on the segment 1 < j < L is 
given by 

L 

Q = Y,c\ Cl . (B.4) 



8=1 



Then we can rewrite (|3j) as 
\(A) A' A >-< 



e lA - l)c\c lt 



(B.5) 



v i=l 



Here we used the projector property (cjcj) 2 = cjcj. Applying the Wick theorem to (B.5) 
we obtain 

X(A) = detT L = det(l + (e' A - l)ff) LxL , (B.6) 
where 1 is the unit L x L matrix and g is the L x L matrix with the matrix elements 



9ij = 9i-j given by (B.3). 



The determinant (B.6) is that of the Toeplitz matrix Tl (its matrix elements (Tl) 



depend only on the difference of the row and column indices i — j). It is said that this 
Toeplitz matrix has the symbol f(e %9 ) with 



f(e l 



l + (e l 



i)e(A* - 



so that the matrix elements are given by 

de 



L)ij 



2tt 



e iq(i-j)f( e V 



(B.7) 



(B.8) 



The formulas (B.6) and (B.3) express the FCS x(A) as the determinant of a given Lx L 
Toeplitz matrix. As we are interested in the limit of large length of a segment, we need 



to find the asymptotics of the Toeplitz determinant (B.6) as L — > oo 



Appendix B.2. The Fisher-Hartwig conjecture 



Let us calculate the asymptotics of the Toeplitz determinant (B.6) as L — > oo using the 



Fisher-Hartwig conjecture [15} 15]. This Appendix is not self-contained. We explicitly 
refer to the notations and the results of Ref. [5J (except for the matrix size denoted L 
in our paper and n in Ref. [9] ) . [§] 



The symbol (B.7) of the Toeplitz matrix has two singularities. In the notation of 



Ref. [5J, one of the singularities must be at 9 = 0, so if one follows the notation of that 
paper, one needs to consider an equivalent problem with the symbol 



it)\ 



for < 9 < 2k F 
for 2k F < 9 < 2vr 



(B.9) 



For a recent similar application of the Fisher-Hartwig conjecture see [THl E] . 
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The two point singularities (Fisher-Hartwig singularities) on the unit circle are 
located at Zq — 1 and Z\ = e 2tkp . These are pure phase discontinuities characterized by 
«o = Oi\ = and j3\ = —f3 = X/2ti (for notations see Ref. [9]). The regular part of the 
symbol f(z) is given by Vjt = for k and by Vo = %2kp^-. Then the theorem 1.1 of 
[H] (due to Ehrhardt [Sj) gives 



X(A) 



(2L sin A; F )" 2K " [G(l + k )G(1 - k )] • 



(B.10) 



Here k = and the result (B.10) is valid for — 7r < A < n in the asymptotic sense as 



L —7- oo. The relative accuracy of (B.10) is o(l) 



The Barnes G function used in (B.10) is defined as [13 



G(l +z) = (2ny/ 2 e-^ E+1 ^ 2 H { 1 + ^) k e- z+ " 2 



k=l 



where 7^ ~ 0.57721 ... is the Euler-Mascheroni constant. We have from (B.ll) 

v 2 \ k 



G(l + z)G(l 



-(7B+l)z- 



n 

fc=i 



e fc . 



(B.ll) 



(B.12) 



To extend the result (B.10) to A outside the interval (— vr,7r), one has to introduce 



the realizations of the symbol 0. In our particular case it amounts to using (B.10) with 
k replaced by 

A 

Kj — 7 

for — 7r + j <\<n+j,j = 0, ±1, ±2, .... In other words, for a given A, one should 



(B.13) 



replace in (B.10) kq — > Kj with the value j minimizing /C. It is important that for A 



being an odd multiple of w there are two such values of j minimizing k -. In these cases 



one should replace (B.10) by the sum over corresponding realizations. In particular, 



according to the theorem 1.13 of jH] we have for A = n 



X(A) 



(B.14) 



£ e 2tkF ^ L {2Lsink F )' 2 ^ [G(l + Kj )G(l - Kj )f 
3=0,1 

Here k .i = ±1/2. We argue in the main text that in fact, the accuracy of this formula 
is 0(1 /L). Moreover, the prescription of adding the most relevant realizations produces 
correct sub leading terms even away from A = tt. 



Appendix B.3. Limits A — > and A = n 

For completeness, we present here the two limiting cases: A — > and A = tt. 



For small A, the leading realization in (B.14) is given by j = (with the term j = 1 
suppressed compared to j = by about l/L 2 ). We have 

X (A « TT) W e ^oL^ln[2Le^ +1 sin fc/ ]-(^) 4 C(3) ? (R15) 

|| A similar improvement of Fisher-Hartwig formula was shown to work numerically for some Toeplitz 
determinants in Ref. [19] . 
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where we have expanded in small z, 
G(l + z)G(l-z) ^e 



-( 7B+ l)*2-iC(3)**H 



We extract from (B.15) 

(Q) 



((Q 2 )) 
((<3 4 >> 



7T 
1 

7T 2 



L = p L, 

In [2Le jE+1 sin k f ] +o(l), 

3 
2^ 



4 C(3) + o(l). 



1,2 



(B.16) 

(B.17) 

(B.18) 

(B.19) 
(B.20) 



((Q 2n+1 »=o(l), n- 

Let us now consider the case A = 7r. Then the dominating parametrizations are 
j = 0,1 (or, equivalently, k = ±1/2). We sum over these parametrizations in (B.14) 
and obtain 

2 cos(k F L) 



X (X = tt) ~ 2 [G (1/2) G (3/2)]' 



V2L sin fcp 



(B.21) 



Appendix C. Details of numerical calculations 



Numerical calculations of the Toeplitz determinant (B.6) 
performed using the superfast algorithm of Ref. 



up to L = 5000, were 
20] implemented in C++ and, 
independently, with the help of Wolfram Mathematica [21] (the two methods agree 
within the error bars, which range from 10 -9 to 10~ 12 , depending on the method). 
For our computations, we took three values of kp (kp/ir = 1/30, 3/17, and 8/17) and 
various values of A (in multiples of 0.l7r). The obtained determinants were then fitted 



according to Eq. (11), and the coefficients Fi(kp, A) were extracted for each pair of the 
parameters kp and A used. This fitting procedure involved splitting the whole range of 
values of L into intervals of an adjustable length and fitting within each interval (using 
Fi(kp,X) and Fi(kp,\ — 2n) as the two fitting parameters). Then the fit parameters 
were extrapolated to L — > oo by using quadratic or cubic polynomials in L^ 1 . This 
fitting procedure allowed us to obtain a very good precision for Fi{kp, A). 

Within the error 



The resulting values of Fi(k F ,X) are presented in Table CI 



bars, Fi(kp,\) is purely imaginary, which, in combination with the reality condition 
Fi(kp, — A) = Fi(kp, A)*, implies that Fi(kp, A) is odd in A. 

We further observe that thus extracted function F 1 (kp, A) can be described by the 



formula (13) to a very high precision. The ratio of the numerically extracted Fi(kp, A) 



to the analytical conjecture (13) is plotted in Fig. CI 



inally, we calc ulate the error e of our expansion at the order L 1 
we show a typical plot of the ratio \e\L 2 



as defined in 



Eq. (|10|). In Fig. 
The 



C2 



Xi(A)| + |xi(A-2tt)| 



plot suggests that this quantity remains finite at L — > oo, which supports our 



conjecture (12). 



In principle, the precision of our numerics should be sufficient for continuing the 
expansion in L -1 (the next term being of the order L~ 2 ). This possibility is also 
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Figure CI. The ratio of the numerically found value of Fi(kp,X) (reported in 



Table CI) to the analytical conjecture (13) 




1000 



5000 



Figure C2. A plot of |e| L 2 (|xi(A)| + |xi(A — 2?r)Q 1 for k F = tt/30 and A/tt = 0.8. 
Here e is the deviation from our conjectured formula, as defined by Eqs. ( 10 ) and ( jll 
with the values of F\{kp, A) obtained from the numerical fit (as reported in Table [CI 
The plotted quantity remains finite in the limit L ~ > oo, which supports our estimate 
p). 
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vy • vy uut: i \ j-j j 


0028957(4)* 


0.6 


51378(3) 

UiUlU 1 VJ \ VJ / 


087217(3) 


0050038(1)* 


0.7 


81586(5) 


138495(3) 


007947(1) 

vy • vy vy i \j jt i i x i 


0.8 


1 21783(9) 

X . — . .1. 1 VJ 1 J 


206731(4) 

vy • vy vy i vy -L i jt ; 


011862(1) 

vy • vy -L j. vy vy j-j i J- / 


0.9 


1 7340(2) 


294348(5) 


016889(2) 


1.0 


2 3785(3) 


403763(7) 

vy • ji vy vy i vy vy i i / 


023166(3) 

vy • vy vy j_ vy vy i vy i 


1.1 


3 1658(5) 


53740(1) 


030832(5) 


1.2 


4 1100(9) 


69770(4) 


04003(1) 

W . W i W W VJ \ X I 


1.3 


5 225(2) 


8870(1) 

vy • vy vy i vy \ _i~ / 


05088(5) 

vy • vy vy vy vy vy \ vy / 


1.4 


6.525(4) 


1.1078(3) 


0.06357(3)* 


1.5 


8.025(9) 


1.362(2) 


0.0782(2)* 


1.6 


9.74(3) 


1.654(1)* 


0.095(4) 


1.7 


11.6(3) 


1.99(2) 


0.12(2) 


1.8 


13(3) 


2.3(1) 




1.9 









Table CI. Numerical values of iF\{kF, X) extracted from the fits. For most of the 
data, the determinants were calculated with the precision 10~ 9 . For several sets of 
data (marked by the asterisks), a higher precision (10 -12 ) was used, with the help of 
Mathematica [21 . Using higher precision results in reducing the error bars for the 
fitting parameters. The dashes in the last lines of the table correspond to the values 
of A, for which F\{kp^ A) could not be reliably determined from our computations. 



connected to the question about the precision of the "sum prescription" (|9j): whether 
this decomposition into a sum (involving all shifts of A by multiples of 27r) is exact to 
all (perturbative) orders in L (see Appendix E) or breaks down at a certain order. We 
leave this interesting mathematical question for future studies. 



Appendix D. Second and third cumulants of the fermionic number on a 
line segment 

In this Appendix, we calculate analytically ((Q 2 )) and ((Q 3 )) up to the orders L~ 2 and 
L _1 , respectively. The Wick theorem expresses those cumulants in terms of the Green 



function (B.l). We further denote the set of lattice sites belonging to the segment 
considered by [L] = {1, 2, . . . , L}. 

For the second cumulant, one gets 

((Q 2 )) = £ 9iiSa (D.1) 



ije[L] 
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where 

9ij = (Qct) = Sij - gji (D.2) 
Using the relation 

oo 

Y 9ij9jk = 9ik, (D.3) 

j=-oo 

one can reexpress 

((Q 2 ))= Y,9QT». (D.4) 

<e[i] 

Collecting together all terms with the same distances between the pair of sites, we arrive 
at 

oo 

((Q 2 )) = 2j2[9,] 2 s(x), (D.5) 

x=l 



where we define the function s(x) = min(x,L). Using the explicit expression (B.3), we 
find 



((Q 2 )> = ^lnL+l[ 7 £ + l + ln(2 sin k F )} 
1 1 cos(2/cpL) | ^, /r _ 3 x 



+ 12^-^(2Limt# + 0(L > (D ' 6) 
(note that there are no terms of the order L^ 1 ). 

Repeating the same procedure for the third cumulant, we find 

((Q 3 )) = ^2 (9ik9ij9jk- 9ij9jk9ik) , (D.7) 

i,j,ke[L] 

which, using relation ( |D.3[ ), can be converted to 

((Q 3 )) = 9i]9jk9u- 9i]9jk9ki- (D.8) 

je[L] i L fce[L] 
3 Mm 3tm 

Again, collecting together all terms with the same relative positions of points, we arrive 
at 

oo 

((Q 3 )) = 6 9x9y9x+y [s(x) + s(y) - s(x + y)\ . (D.9) 

x,y=l 

Note that there is no contribution to this sum from small x and y, for which x + y < L. 
At large L, the main contribution to this sum is determined by the two boundary pieces: 
(y = 0, x > L) and (x = 0, y > L). The easiest technique to extract these contributions 
is summation "by parts" . If we re-factorize 

,, r n s(x) + s(y) — s(x + y) _ 

((Q )) = 6 Yl [xy9xg y gx + y] ■ — — — x — - — - , (d.io) 

x+y>L 



Quantum fluctuations of one- dimensional free fermions 



15 



then we can associate the leading contribution to the jumps of the function [s(x) + s(y) — 
s(x + y)]/{xy) at the boundary. The contribution from the two pieces of the boundary 
are equal, and, after a simple calculation, one arrives at 

3 ^-^ sin(2k F y) 

73 



«Q 3 » 



vr- ^ x(x + y) 

y>0 

Af>n(2^).| L % 

y=o JL 



2tt 3 L 



cot ki 



(D.ll] 



to the leading order in L (in the calculation, we neglected the terms oscillating in x in 
QxgyQx+y, as they produce contributions of higher orders). 

Appendix E. Charge cumulants from the generalized Fisher— Hartwig 
formula 



In this Appendix, we try to push our conjecture about the precision of the generalized 
Fisher-Hartwig formula even further than in the main body of the paper. Namely, we 
assume that the generalized Fisher-Hartwig "sum prescription" is exact to all orders of 



L 1 , if all possible realizations (B.13) are included: 

+0O 

X(A)= Yl X*(A-2ttj), 



(E.l) 



where 



X*{\) = exp 



■\ kp T 



X 2 L 
2^ 



oo 

£ 

m=0 



F m (k F ,X)L~ 



(E.2) 



(we assume that the sum in the exponent contains only L m , but no logarithms). The 



series (E.l) may be, in general, divergent at a fixed L and should be understood as an 



asymptotic expansion at L — > oo. Such a conjecture may, in principle, be verified by 
comparing with the charge cumulants (|4]). In particular, for the first several cumulants, 
we find from Eqs. (E.l) and (E.2j^: 



(Q) 
((Q 2 )) 

((Q 3 )) 



-k F L . 



^2LsinM-4 r5 2 ^ L l 



7T 

— ln(e^ Tj-2 (2Lsin A;^) 2 

- F'(L- X - F^L- 2 + o(L- 2 ) , 
6 sin(2k F L) 



+ 



ln(e 7E 2L sin k F ) 



vr 3 (2Lsmk F ) 2 
+ iFf ) L- 1 + tF^L- 2 + o(L- 2 ) 



((<3 4 >> = - 



3 24 cos(2k F L) . 2 

2^ C(3) + ^ (2LsinM 2 N^2LsmM] 



+ F^L- 1 + F^L- 2 + o(L 



(4) r -2 



-2\ 



(E.3) 



(E.4) 



(E.5) 



(E.6) 



% Here Fm' means d*F m 



A=0 
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In other words, the coefficients in the expansions of the cumulants ((Q n )) in L _1 may be 
related to the expansions of F m (kp, A) in A (at A = 0). These relations may, in principle, 
be verified against direct calculations of the cumulants using the Wick theorem (along 
the lines of Appendix D). 

In the present work, we have only calculated analytically the cumulant ((Q 2 )) to 
the order L~ 2 and the cumulant ((Q 3 )) to the order L^ 1 [see Eqs. (D.6) and (D.ll), 
respectively]. This implies = 0, F% = -1/(12tt 2 ), and F x (3) = —(3i/2ir 3 L) cot kp [in 
agreement with our conjecture (13)]. Furthermore, numerical studies of the third and 



fourth cumulants indicate that ff = F^ = and F?> ± 0. In other words, if our 
conjecture in Eqs. (E.l) and (E.2) is correct, then F 2 (kp,X) has the form 

F 2 (k F ,X) = -^ + 0(X 4 ) 



,(4) 



(4) 



247T 2 



(E.7) 



in its expansion around A = 0. It seems likely that further analytical progress is possible 
in these questions, and we leave them for future studies. 
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